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The cross section value of doubly charmonium production in e + e~ annihilations 
have been estimated at interaction energy ys = 10.6 GeV. The quarkonium wave 

O ' 

| function shape, as well as the non-zero value of charm quark mass, have been taken 



< 



(N 



X 



into account. 

The recent experimental results of the BELLE and BABAR Collaborations on double 
charmonium production in e + e~ annihilations demand an essential revision of the calcu- 



lation techniques based on QCD factorization theorem. Indeed, the cross section value of 
charmonium pair production in e + e~ annihilations estimated within standard calculation 
; ■ technique underestimates the experimental data by an order of magnitude 

a a. The 

detailed analysis has shown that one of the reasons of such underestimation is large fixed 
virtualities of the intermediate quark and gluon, q 2 ~ ^- (Q 2 is a virtuality of the initial 
photon, see Fig. Q]), which occur in the standard calculations, where relative momenta of 
valence quarks of charmonium are not taken into account in the hard part of amplitude (the 



so-called (^-approximation). 

Taking into account the relative motion of valence quarks allows to decrease contradictions 
between theory and experiment. Calculations of thepair quarkonium production in e + e~ 
annihilations, obtained within a light cone formalism jj, O, are in a qualitative agreement 
with the experimental data of the BELLE and BABAR Collaborations. In these works, the 
5-shaped wave function of quarkonium was replaced by a function (j)(x), which is "spread" 
on x, where x is a momentum fraction of the quarkonium carried by a valence quark in the 
infinite momentum frame. Thus, the longitudinal "internal motion" of the valence quarks 
is taken into account (see also pi]). This spreading decreases the effective virtuality of the 
intermediate quark and gluon, and, therefore, increases the predicted cross section value. 

An analysis of doubly charmonium production within quark-hadron duality also leads to 
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Figure 1: The leading order diagrams for the 
process e + e~ —> J/ip + cc. Q is momentum 
of the photon and q is momentum of the in- 
termediate quark or gluon. The average of 
q 2 in the total phase space is about several 
m 2 . If both cc-pairs are fused into charmonia 
(e~e + — ► J/ip + charmonium), then q 2 ~ Q 2 /4 
within the ^-approximation. 



Figure 2: The shape of charmonium wave func- 
tion used in J] within the lightcone formalism 
(leading twist). 



the cross section value increase [8], and allows to describe the data. 

Recently several works have appeared, in which the internal motion effect has been esti- 
mated in the framework of NRQCD [si, Q. The received v 2 corrections increase the cross 



section by a factor a 



cone approximation [J, [5. 



Dout 1.5 — 2.5. This value is less than value ~ 10 obtained within light- 

rn 



where the critical analysis of lightcone 



]. Also there is a study 

formalism has been done for the double charmonium production in e + e~ annihilations. The 
predictions obtained in this work are close to the NRQCD estimations 

In this work we study the dependence of internal motion effect from the wave function 



shape, as wel 
described in 



as from the c quark mass value. For this purpose we use a simple model 



12j instead of rigorous lightcone formalism or NRQCD. Within the model 



under consideration, an amplitude of J/iprj c production in e + e annihilations is performed 
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as follows: 



Aj/^nc = J J dx dy (f)j/ip(x)(j) Vc (y)T cSc5} (1) 

where <frj/^(x) is a wave function of J/ip, x is p + fraction of c quark inside J/ip, (pr] c (y) 
is a wave function of r], y is p + fraction of c quark inside r] c , and T c5c c is a hard part of 
amplitude, which describes the process of four charm quarks production. The cc-pairs in 
T C ccc have quantum numbers corresponded to J/ip and rj c states. The transverse motion of 
quarks inside charmonium is neglected. 

For <f)j/ip and <f) Vc functions one choose the parametrization used in work |4| within the 
lightcone formalism: 

n 1-v 2 



<Plightcone{x,V 2 ) = c{v 2 )<f) a {x) 



x(l — x) 



(2) 



1 -4x(l -x){\ -v 2 ) 

where c{y 2 ) is a normalization coefficient, <p a (x) = 6x(l — x) is an asymptotic form of wave 
function and v is a typical velocity of quark-antiquark pair inside the charmonuim. In accord 
with the potential model for J/ip and r] c mesons v 2 ^0.23 H. 

The results of these estimations have been compared with the predictions of the 5- 
approximation v 2 ) = 5(x — 1/2)) and with the prediction for the case of "tailless" 
wave function: 

2 (2.5 0.3<x<0.7 

<f>tailless{X,V ) = < (3) 

[ x < 0.3, x > 0.7 

The cross-section estimations within the model (DQ) for the wave functions ((2j) and ([3D 
have been done numerically. The standard expression for 5-approximation is 

(1 - AM 2 /sfl 2 
a , (4) 

where s is the e + e~ interaction energy squared, and M 2 is a charmonium mass. 

The ratio between our calculation results and the ^-approximation are shown in Fig. [3] 
for different values of charm quark masses and wave function shapes. The ^-approximation 
predictions are given for m c = 1.5 GeV. 

The toy model under discussion do not allows to shed light on all aspects of double 
charmonium production; nevertheless the two main features of such process are clearly seen: 
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Figure 3: The ratio between the calculation results and the ^-approximation for different values of 
charm quark masses and wave function shapes as a function of interaction energy: upper solid curve 
is for m c = 1.25 GeV and the lightcone wave function (|2j); lower solid curve is for m c = 1.5 GeV 
and the lightcone wave function j2j); upper dashed curve is for m c = 1.25 GeV and the "tailless" 
wave function ([3]); lower dashed curve is for m c = 1.5 GeV and the "tailless" wave function (J3j) . The 
^-approximation predictions are given for m c = 1.5 GeV. 



1. The cross section value is very sensitive to the c quark mass value at interaction 
energies below 15 GeV. 

2. The cross section value depends strongly on the wave function behavior on the tails 
(at x < 0.3 and x > 0.7). 

As one can see from Fig. [3J the uncertainties caused by c quark mass value are large at 
energy 10.6 GeV, where Belle and BaBar data on doubly charmonium production have been 
obtained. The cross section value for wave function ((2j) is increased by factor 1.9 for m c = 
1.5 GeV and by factor 3 for m c = 1.25 GeV in comparison to the 5 approximation. For the 
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"tailless" wave function (J3]) the cross section enhancement is about 20% for m c = 1.5 GeV. 
For m c = 1.25 GeV the cross section value is increased by factor 1.6. Unfortunately the 
tail behavior is not known with sufficient accuracy. Also it is not clear which c quark mass 
should be used to obtain the cross section. For the (^-approximation, it should be better to 
use m c = M/2 to obtain the correct meson mass. The situation is not so obvious if one takes 
into account the quark internal motion, because this motion increases the invariant mass 
of quark-aniquark system. This a reason to choose the c quark mass value slightly smaller 
than M/2. It is worth to note, that for the "massless" c quarks the model under discussion 
leads to the cross-section increase by factor 10, and, therefore, reproduce approximately the 
predictions of lightcone formalism. 

Contrary to the double charmonium production in e + e~ annihilation, the internal motion 
can be neglected for the process of associative charmonium production e + e~ — > J/tpcc. The 
cross section estimations for this process have been done within the same model, as for the 
double charmonium production: 



Our estimations show that the internal motion changes the cross section value by 5-10% in 
comparison to (^-approximation [lj|. These results are in agreement with predictions in flo| . 
The dependence on c quark mass is more essential for this process. 
To summarize: 

• Within the model under discussion the cross section value for the process e + e" — > 
J/iprjc at 10.6 GeV can be increased by factor 1.2-3 depending on the c quark mass 
value and the wave function shape. 

• In the frame of this model the dependence on the wave function shape can be neglected 
for the process e + e~ — > J/ipcc. 

• The additional theoretical studies of charmonium production processes in e + e~ anni- 
hilations are needed. 
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